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Introduction 


A Spacecraft Control Laboratory Experiment (SCOLE) has been proposed by 
Taylor and Balakrishnant to evaluate control laws for flexible spacecraft. 

The SCOLE configuration consists of a large antenna attached to the space 
shuttle orbiter by a flexible beam. Taylor, et. al., discusses the details of 
the experiment and offers a mathematical model of the system dynamics. The 
dynamics are described by three distributed parameter beam equations with rigid 
body attachments at both ends. The boundary conditions at the beam ends 
contain the forces and moments, both inertial and applied, of the rigid shuttle 
and reflector bodies. Non-linear kinematics couples the otherwise uncoupled 
beam equations. 

This paper presents the linearized (dynamics of the SCOLE configuration. 

The equations are developed by the method of Lagrange and are assembled into 
Matrix Second Order Form. When appropriate, the nomenclature, sign conventions 
and coordinates used by Taylor, et. al., are adopted for convenience. 


t Taylor, L. W. , Jr. and Balakrishnan, A. V., "A Mathematical Problem and a 
Spacecraft Control Laboratory Experiment (SCOLE) Used to Evaluate Control 
Laws for Flexible Spacecraft... NASA/IEEE Design Challenge" (Unpublished). 
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Nomenclature 


Coordinate Systems 


1 ' 

A A A T 

[i, k] 1 

Inertial Frame 

2 " 

^b* J b’ y T 

Shuttle body-fixed Frame 

3 ~ 

[5 e > 3 e , k e ] T 

Elastic body (beam) Frame 

4 ‘ 

A A A T 

[ipj k r ] 

Reflector body-fixed Frame 


Transformations 

Tjj - Inertial to shuttle body, i.e., ^ = T^Ij (T^ in [1]) 

Tfar" Reflector to shuttle-body, i.e., ^ = T br I 2 (T^Tj ln [1]) 

Symbol s 

$ - total displacement of a point on beam 

C 

- total displacement of reflector c.g. 

? r - vector from shuttle body frame to reflector frame 

r b - vector from inertial frame to shuttle body frame 

s - distance along beam in z g direction w.r.t. Ig 

a b - vector from shuttle body frame to beam frame 

a r - vector from reflector body frame to beam attach point 

f - proof mass force vector 

F r - force vector at reflector c.g. 

r g - vector from beam frame to any point on the beam 

T, - control moment applied to shuttle body 
s 

T r - control moment applied to reflector body 

Tp - disturbance torque applied to shuttle body 

J - product of inertia of shuttle body (w.r.t. c.g.) in 
s xz the xz-plane 
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- moment of inertia of shuttle body (w.r.t. c.g.) about x,y,z axes resp. 

- moment of inertia of reflector body (w.r.t. c.g.) about x,y,z axes resp. 

- shuttle mass 

- reflector mass 

- beam mass 

- beam mass density . cross-sectional area 

- pitch angle, rotation about y axis 

- roll angle, rotation about x axis 

- yaw angle, rotation about z axis 

- deflection of beam in y-z plane 

- deflection of beam in x-z plane 

- angular deflection of beam about z-axis 

- length of reflector mast, beam 

- bending rigidity of beam in y-z plane 

- bending rigidity of beam in x-z plane 

- torsional rigidity of beam about z-axis 

- Mass Matrix 

- Stiffness Matrix 

- Disturbance Input Matrix 

- Control Input Matrix 

- acceleration, attitude and angular velocity sensor noise 
respectively. 

- plant noise 

- i tJl generalized coordinate 

- measurement vector 

- beam mass density * beam polar moment of inertia 

- beam polar moment of inertia • beam cross-sectional area 
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Kinematics* 


The total displacement of the shuttle body at the c.g. is defined 
by. 


(la) 



and the total displacement for an arbitrary point on the beam and at the 
reflector c.g. are defined respectively by 

(lb) fi e (s, t) = + r e = l]r b + lja b + lju(s, t) 

= X l^ + T b ( ^b + t ®** 

where u(s, t) T = [-u 0 (s, t), U (j> (s, t), -s]*** 

(lc) * r = r b + r r = ijr^ + ija^ + lJu(L, t) - ija^ 

= l{[r b + Tj(a b + u(L, t) - tJ^)]** 


T 

For small angles 8 and ip the transformation matrix T b is 
defined by. 




1 

rp 

-e 


1 


0 

- 4 > 

1 


and the transformation matrix T br by, 


*see figure 2. 

**l2 and are coincident 

***The mast is treated as an Euler-Bernoull i rod and u , u and 
u, are assumed small. * 6 
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'br 


1 

-%< L ' 

u'U, t f 

V L * t> 

1 

-uJ(L, t) 

-Ug(L, t) 

uJ(L, t) 

3U(.)(L» t) 
Ss 

1 

h_ 

\(L, t) = 



For both and T^, the NASA-Standard Euler Angles for rigid body 
rotation are assumed. 

The velocity equations: follow directly from equations la,, lb and 
lc. 


(2a) 


? = I T f 
r b ‘l r b 


(2b) i5 e (s, t) = l{[r^ + + iju(s, t) + lju(s, t)] 

(2c) fl r = ifo + + u(L, t) - tJ^) + Tj(u(L, t) - tJ^)] 

Making the proper substitution equations 2a, 2b and 2c become. 


(3a) 


VI 


(3b) ff e (s, t) = I. 


“l 

x 

y 

z 


x-^-se-Ug-^u^ 

y4u 0 +s$-iHi e +iy 
• • • • • 
L z+eu e +$u^+eu e +<i>u ( ^ 


where u^.j = u (.)( s » *) 
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(3c) 9 « W e (L, t) + i] 


(<,„ i)i+ eu^*u i|j +6u;) a r ^(u 4i -6u | ;+*- e ^ ) a l 

• • • • • • •• 
-(’i'-M'U'+u + 4 < )a„ +(^u +<,u!+^u + 4 >u!)a 


e g> ' e r 


4> Y ♦ v > 4> r 


. (e -' > % +u e‘% )a r -( 6 V‘» + eu | +u’)a r 
x y 


where u^.j = U (.)( L * t). 


By assuming that products of small terms are neglible, equations 
3a, 3b and 3c simplify to, 


(4a) 


-v t T 
r b = ! 1 


X 

y 

z 


(4b) ft p (s, t)=l{ 


X-S0-U 


e 


y+s<j>+u 


♦ 


(4c) 


ft =1 

w r X 1 


x-Le-u 0 +(u ij; +^)a r 

y +L i + %-(i%)a r 

z+(e+uj)a -(i+u^ a r y _ 


“(•) " U (-) (L ’ 


Kinetic Energy, Potential Energy and Generalized Forces 
In the method of Lagrange, the equations of motion are formulated 
from the scalar quantities of kinetic energy, potential energy and 
work. Using the kinematics of the previous section these three 
quantities are easily written; 



Kinetic Energy 


T = l/ 2 {M_(x 2 +y 2 +z 2 )+J e i 2 +J e e 2 4 J c i 2 + 2 J C U 


x y z 

( 5 a) +M (w 2 +W 2 +W 2 )+J (i+u*(L, t )) 2 
x y z x 9 

+^ r (e+u'(U t)) 2 +J (i+u (L, t )) 2 
y 'z 9 


xz 


L 

r < 

0 


x y z 


+ / o(s)(W 2 +W 2 +W 2 )ds + / <(s)(fHi (s,t)) 2 ds} 
e ~ e.. e_ 0 v 


Potential Energy 


( 5 b) U = 1/2 / L {EI e u" 2 (s, t)+El u" 2 (s, t)iGJu; 2 (s, t)}ds 

o 99 


Work 


( 5 c) 


4 W = iQ.lq. * <F_ +f, +f 2 >V< F r +f 1 +f 7 
-! 1 1 x x x * y y y 


+ (0)5z + ( T s e +T re +T D e " f i x s r f 2 x s 2' F r x L)66 

+(T S +T r +T 0 +F L+f. s +f s ? )5* 

♦ * ♦ y y y 

+( T S +T +T 0 -F a +F a )«* 

* ■ ♦ * x y y x 

+T r Su^L.tMj 6u e (s is t)-f 2 6u 0 (s 2 ,t)-F r <5u 0 (L,t) 


+T 6u ' (L,t)+f, <Su (s,,t)+f 9 Su (s ? ,t)+F 6u (L,t) 
4 > i y 4 > 1 d y <p t r y 

+ T su (L,t)-F a Su (l.O+F a «u (L.t) 

\p 9 x y 9 y x . 


Defining the following Ritz approximations for u^, u 0 and u f 
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0 

u e ( s » t) = Jv. (s) n . (t) 5 

ie ; 

Ms. t) = £v. (s)n. (t) 

v 1 * U 

n 

Up(s. t) = £ (s)n (t) 

and making the appropriate substitutions the expressions for the kinetic 
and potential energy and work. can be written, 


T = l/2M s (x 2 +y 2 +z 2 )+l/2(J e 2 +J $ 2 +J * 2 ) +J 

s y s x s z s xz 

+1/2M (i-Le-Jv. (L)n, (t)+k. +R. (L)rj. (tU ) 2 
i V ’e r y i \ \ r y 

+1/2M (y+L$+£v. (L)n- (t)-rf-a -]Tv. (L)n. (t)a ) 2 
i * 4> r x i \ \ r x 

(6a) +l/2M r (z+ea -*a +£v! (L)n, (t)a -£v! (L)n- (t)a r ) 2 

x y i 6 'e r x i \ > r y 

+1/2J (e+Jv! (L)n, (t)) 2 +l/2J (i+Jv. (L)n. (t)) 2 

y i e e r z i ^ > 

+1/2J U+£v! (L)n. (t)) 2 +l/2 J L o(s)z 2 ds 

x i 1 4> 4 > o 

+1/2 / o(s)(x-se-Tv. (s )t» - (t)) 2 ds 
o i ^ ’e 



+GJ(Iv'. (s) n , (t)) 2 )ds 
i \ \ 


(6b) 



= (F +f, +fp )6X+(F +f, +f 2 )6y 
“ “ “ y y y 


x x x 


+(0)6z+(T +T +Tq Sp-F L ) 6 0 

e e e x x x 

+( \ + VV F r„ L+ \VV2 )S * 

4> 4> y y y 

+(T +T +Tpi -F 3 +F 3 ) 6\i> 

' s r D r r r r y 

x y y x 

+I(T r v! (L)-f. v. (s,j-f 2 v. (L)-F v. (L»6n i (t) 
i r e * e Sc x r x ’e *e 

+ I ( Vi {L)+f l . v i ^ s l )+ ^2 v i , ^ S 2 )+F r v i (L))5n i {t ^ 
■» 4> <t> y <t> y <*> y $ 4> 

+7{T_ V. (L)-F a V {L)+F a V <L))«n. (t) 

i t * x y * y x * > 
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Lagranqe*s Equation 
Lagranges Equation is given by, 

* 7 - * __ 

d /3T \ ST , 5U _ « 

- aq f - V 

where q. are generalized coordinates, 


Defining the generalized coordinates as, 


[q^» i=l **p]-[x,y,z,$,e,^,n^ (i=l-n ), (i=l-»n e ), r\. (i=l-*n. )] 

<J> * B |J; ■ ™ 

where p=6+n+n a +n, 
r $ 0 


and performing the required operations on equations 6a and 6b and using 

\ 

5c results in the equations of motion shown in Appendix I. Assembled in 
Matrix Second-Order Form, the equations can be written, 

Mq + Kq = Bu + Dw* 



% 

% 


0 0 


V 

where M = 



, K = 


D = 




V 




0 


u T • [F , F , f, , 
x y x 

w = [T d » T q , Tp ] 

4> Q 



A numerical example is presented in Appendix 3. 



M 


RR 


M +M +m. 
srb 


0 

0 

0 

m h L 

-V--T 


M +M +m K 

STD 


0 

0 




m b L 


m. L " 14 ' 

-V-f Vr. . 


r 2 


m b L 


fMJ. 


2 r 

a 


M+M+m, -Mb 
s r b r r. 


-M a 
r r. 


0 


M a 
r r 


-M a 
r r. 


“ M r a r a r J s “ M r La r 
r r x r y s xz r r x 


M a 
r r. 


-M -a a 
r r r 
x f y 


'W 


Vr 

L y 


'Vr 


,0 


J s - M r La r 
XZ X 


- M r L V 


J, = 


Jo = 


Jo = 


J c +J„ +M L 2 +H a 2 +m.L 2 /3 
s x r x r r r y b 

+J r +M r l 2 +M r a r +m b L2/3 

y y x 

J +J +M (a 2 +a 2 )+J. 
s r r r r b 
z z y x z 
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$.y 

\.y 

4 
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•*1 a r 
$.y x 

■ XJ .A 

'*1 

$.x 

-X 2 

*.x 


M RE T 


0 

‘ 0 

X* 

X 1 

e.x 

x 2 

e.x 

- X 1 8 r 
e.y r x 

-*2 “r 
e.y r x 

V 

4 

e.x 

-*1 
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* X 2 
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x i 

- X 1 

4<,x 

-X 2 
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• 

0 

-“1 
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S-.x 

• 
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*.e 

-b 12 

*.e 

" 6 n 

4>.4> 

~*12 

4>.4> 

' 2 >../ 2 >. 

,22 ..; p22 . 

~®21 

**.e 

_ 0 l • • 

b 22 

*.B 

. • • 

- 6 21 
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*®22 

♦ .4- 

•*ii 

e.* 

"®12 

e.4 

Xn e.x +0ll e 

• • 

. ■ • • 

X.x 12 . 

-Xj Vj (L) 
e.y ** 

•X-i v, (L) ••• 
e.y % 

”®21 

e,$ 

-e 22 

e,$ 

X.*' 21 . 

Xx*' 22 e 

"X? V 1 
e.y * 

-x 2 v 2 (L) 
e.y 41 

" 6 11 

-6 12 

>('•4 

*Xj vj (L) 
«.y e 

-Xi V? 

4-.y e 

e ll 

E 12 
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Y i v * M r v i ( L ) + / (s)ds 

’(*) r V) o '[•) 
yj) = M a J) (L) 

’(•k-o r r (-o v> 

P M(.j 4 


Note: For compactness in 

notation note that J has 

♦ .e,* 

been used in lieu of 

r x,y,2 




8„ * M a v' (L)v'. (L)a 

”(•)(••) r r y ’(•) >(■■) r 


a, ■ (J r +M a 2 +M_a 2 )v. (L)+ J ^(slv, (s)ds 

’(•) (•) r r y r V ’(•) o ’(•) 

«j = M r La r (l) 

’(•)(••) r (••) ’(•) 

T «(.)(-)‘ ^(••^(•) a) ^(.) (L N-) v n a ) )v i(.) (L)+ V (5 \., (s \.) (s)ds 

5 ’c )(••)* ' V '(-A.) )v V) (L)+ )( l) 

A 2 o L 

e ii = M r a r v i ^ L ^ v i ( L ) +J r v i ( L ) v i (D+M v. 0-)a* v -i 0-) + / ( s ) v i ( s ) d s 

J y \ J * 2 > J t|; r \ r X 0 > 


d r = 


0. 

0. 

0. 

1. 

0. 

0. 


0. 

0. 

0. 

0. 

1. 

0. 


0. 

0. 

0. 

0. 

0. 

1. 
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II 


CO 


“ location of proof mass actuators 



18 


EE 


k ll k 12 

<J> 4> 

n a 

k 21 k 22 

$ 4) 

• • 

• • 

0 0 • • • 

0 0 • • • 


0 0 

o o 


““e k]2 e " 

U It 

e “e 


o 

o 


0 

0 


0 

0 


0 

0 


0 0 

b o 

• * 

» * 

• m 

o 6 

o o 


k k* 

K ll t( K 12 

ip 

k' i/ 7 .. 

21 . k 22 

*1‘ ip 


J 


where 


A 1 5 f (s>vf ( 

J (‘) o < ' '(.) J(.) 


s)ds 


note: for the $ terms, 6J is required in place of EL 

Measurement Equations 

Per Taylor et. al ., both the shuttle body and reflector are instru 
mented to obtain inertial attitude, angular velocity, and acceleration 
measurements for all axes. Additionally, two axis accelerometers are 
located at intervals of 10 feet along the mast. 

The corresponding measurement vector is defined by 




= [<J> S » 6 S » \p s » 4y» 0j,, » Z = 2 


0! A 


= [x s , y s , z s , x r , y r , z r , x b , y b ; i - 1-*12] 


where 
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and (•!» (•) » and (Ok denote shuttle body, reflector body and 
s r d ^ 

location of the i** 1 beam sensor respectively. 

Manipulating .the kinematical equations previously defined yields 

the following measurement equation. 


4 



where. 


z a = 

R „ = P A- z = z - 

“ T=! t w V A V „ V J 



M, K and D as previously defined 
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Line of Sight Calculation 



Figure 3 - Line-of-Sight Definitions 

The line-of-sight (LOS) is defined by a unit vector (\ os ) from the feed 
which is reflected at the center of the reflector.' The expression for 
\oS is eaS1 ^ derived ( se e Appendix 2) and is given in inertial 
coordinates by. 


10S 


(VV-2[(\-R f )-R a ]R a 

|(^ r 4 f )^2[(R r -R f )-R a ]R a | 


where R„ 


= 1 , 


/° 1 
o’ 

1 « J 

\ = I T T T T T 
I *1 T b T bd 

f n 

1 ? 

\ l i 


1 


and R r , Rf are defined in 


the above figure. 
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The LOS error is affected by both the pointing error of the 
Shuttle body and the misalignment of the reflector due to .bending of 
the mast. The LOS error is given by. 


LOS 


error 



where Dy is the taVget direction. For small LOS error angles. 


LOS 


error 




APPENDIX 1 


Translational Equations 


(7a) 


(H s +« r -Hn b )x-(« r L4m b L/2)e-M r Jv d (L)nj (t) 

J 0 6 

L 

- / o(s)jv. (s)n.- (t)ds+M a ^+M a Tv. 

o 3 J e J e y yj 


= F r + f a + f 2 
x x 


X 



(t) 




Rotational Equations 

(J s +J r +M r l2+M r a2 r + m b L 2 /3H+(m b L/2+M r L)y 


-M r a r z-M r a r a 

«y y 




+M r Llv. (L)n . (t)t(J -MJ_a r )i-M La r Iv. (L)ri. (t) 
r j s xz r r x r r xj ^ 


(7d) 
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= VVYV'V^V 2 


( 7 e) (J +J +HL Z +Ha 2 +m b L 2 / 3 )e-(M r L+m b L/ 2 )x 
y y x 






( 7 f) (J +J r +J h +M a 2 +M a 2 )J+M_a x-M a Le 
2 2 2 y x y y 



L 

-M a y+(J c * M r L? r ^‘ M r a r £ v i (*-)n* (*)+/ K ( S )I v. (s)n. (t)ds 
r r x s xz r r x r r xj b J <p 0 j b b 


= b +T r +T D ~ F r a r +F r a r 
b iji r x y y x 


4 

Beam Bending and Torsion Equations 
L 

(7g) (M v. (L)+ / o(s)v. (s)ds)y-M v- (L)a * 
r % o $ b x 

-Ha v' (L) 2 -M a a v' (L)e 

y b y x b 

l 

+{M v. (Dltlli »: (L)+J r v! (l)+ / ots)sv. (s)ds)i 

<j> y b x o V 
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= 1 v! (L)+f. v. (sj+fp v. (s 2 )+F v. (L) 

♦ 4> y 4> y * y ♦ 


L 

(7h) -(M v. (L)+ / o(s)v. (s)ds)x+M v! (L)a z 

e o 0 e x 

-M v. (L)a i-M v! (L)a a * 

e y e x y 

+(M v. (L)L+M v ! (L)a^ +J v! (L)+ / o(s)v, (s)sds)e 

e e x r y e o e 

+M r v. (L)J\> • (L)n . (t)-M v. (L)£v. (L)n. (t)a 

r ’e 3 J e J e r ’e j j m > r y 

+M r v! (L)Ivj (L)n, (t)aj -M v! (L)£v! (L)ri . (t)a a 

r ’e j J e J e r x r ,j J <}> °<j» r x r y 

L .. 

+J r v! (L)[v! (Ljn, (t)+ J o(s)v. (s )£v. (sjn,- (t)ds 

ye j u e J e o e j J e °e 

L 

+ / (t))*? (s)ds 

o j J e J e e 


= T r. v i. (L) - f l. v i. ( ‘l>-Vl.<* 2 ) - F r v i. (L) 

x 


e ’e ‘x 'e “ “x ' e ~ r .. e 
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(7i) H a v. (L)x-M a lv. (l)e-K v, (L]a, y 
rr y't r r y % r ’♦ r x 

-H-v, (L)a U-M a v. (L)£v. (Lh, (t) 
r x y > j J e J e 

+(M r a^ (L)+M r a^ v, (L)+J r v. (L>»-/ k(s)v- (s)ds)>. 
y <l> x tp 2 v ° \p 

+M r a r v i (L ^ v j tL)n j +J r v i (L)n i + / “(sK (s )R. (s)J. (t)ds 

y V* 3 z '* j •> > o j 0 ^ ^ 


~ M r a r v i ( t ) +M r v i t L ) a r l v i ^ L ^ n i (*) 

x > j J 4> • \ r xj 

L 

+’■/ El (|vl (s)d (t))vi (s)ds 
o v j J \i/ P 

- T r v. (L)-F a v. (L)+F a v, (L) 

* v r x y V r y r x > 



APPENDIX 2 



/ \ 

I \ 


From the above figure, 

(D = p+q and ^ L0S * 

It is easily seen that, 

ti -+C(R r -R f )-R a ]R a 


and from (1), 

p= (* r -fy-5 


Then p-q = (R r ~R f ) - 2q 

= (R r -R f ) - 2[(R r -R f )-R a ]R a 
and _> 

j . (yV-2C(R r -R f )-R a ]it, 
105 ' KV*f)- 2 n« r -s f )-s a ]* a i 

—V 

where R 0 has unit length. 
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APPENDIX 3 

In the following example, the beam was modeled with four finite 
elements using cubic polynomial shape functions for bending and linear 
shape functions for torsion. 


the 

mass 

matrix < 

26 by 26 

) 




1 

2 

3 

4 

5 

6 

1 

6# 40e + 03 

0 • 

0 . 

0. 

-2.42e+03 

4«04e*02 

2 

0. 

6* 40e*03 

0 . 

2.42e+03 

0. 

2.33e+02 

3 

0 . 

0. 

6, 40e*0 3 

-4. 04e+02 

-2.33e*02 

0. 

4 

0 • 

2.42e+03 

-4*04e+02 

l«20e+06 

7. 58e+03 

-1.15e+05 

5 

-2.42e+03 

0 . 

-2.33e+02 
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